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Abstract. We compare magnetic field extrapolations from a photospheric magnetogram with the observationally 
inferred structure of magnetic loops in a newly developed active region. This is the first time that the reconstructed 
3D-topology of the magnetic field is available to test the extrapolations. We compare the observations with 
potential fields, linear force-free fields and non-linear force-free fields. This comparison reveals that a potential 
field extrapolation is not suitable for a reconstruction of the magnetic field in this young, developing active region. 
The inclusion of field-line-parallel electric currents, the so called force-free approach, gives much better results. 
Furthermore, a non-linear force-free computation reproduces the observations better than the linear force-free 
approximation, although no free parameters are available in the former case. 
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1. Introduction 

Due to the low plasma f3 the magnetic field is the 
dominating quantity in the low solar corona. Thus, 
the 3-D magnetic field structure is of basic importance 
for physical processes in the solar atmosphere, such as 
flares, coronal mass ejections and X-ray jets. Direct 
observations of chromospheric and coronal magnetic 
fields arc difficult, but significant p rogress has been 



made within the last fe w yea rs. 
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Lagg et all (|2004h ; iLin. Kuhn. fc Coulter! feOoA ). Here 
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we compare magnetic loops reconstructed from magnetic 
field measurements bv lSolanki et alj (|2003[ ) with magnetic 
fields extrapolated from a photospheric magnetogram. 
The measured fields allow a much more sensitive test of 
extrapolations than other observations. 



Vacuum Tower Telescope (VTT). The spatial resolution of 
the data was limited by seeing to 1.5". 

The photospheric magnetic vector map was obtained 
by applying th e inversion co d e SP INOR to the Si I 
Stokes profiles (jFrutiger et all l2000h . The He I triplet 
provided the maps of the chromospheric vector mag- 
netic field. Th is triplet, wh i ch ha s a complex non-LTE 
line f ormation ( Avrett et all 19941 ) but is nearly optically 



thin ( Giovanelli fc HaliTl977 ), was analysed by applying 
the U nno-Rachkowsky solution ( Unnol . 19561 : [RachkowskvL 
1967) to describe the individual Zecman components of 



19671 ) to describe the individual Zeeman components of 
each member of the triplet, together with a simple im- 
plementation of the Hanle effect based on recent devel- 
opments. These have convincingly demonstrated that the 
Hanle effect in forward scattering close to the solar disk 
cente r creates measurable linear polarization in spectral 
lines (iTruiillo Bueno et all [20021 ). 



2. Measurements of magnetic fields in the upper 
chromosphere 

The inference of the magnetic vector is based on an in- 
version technique applied to spectropolarimetric data of 
the photospheric Si I line at 1082.7 nm and the chromo- 
spheric He 1 1083 nm triplet. The data were recorded with 
the Tenerife Infrared Polarimcter mounted on the German 



Send offprint requests to: 
wiegelmannOlinmpi .mpg . de 



T. Wiegelmann, e-mail: 



The inclination and azimuthal angle from the chro- 
mospheric magnetic field map was used to trace magnetic 
field lines. We identified field lines as magnetic loops if the 
following criteria were fulfilled: the magnetic field strength 
must decrease with height, the inclination and azimuthal 
angles must not vary strongly from one pixel to the other 
and the height of the two footpoints must be similar. 
For a more detailed description o f the observations and 
the analysis techniq ue we refer to lLagg et al.l (|2004h and 
Solanki et all (12005 
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We are aware of the assumptions and approximations 
which enter the magnetic field reconstruction from the 
polarimetric observations: The Milne-Eddington approach 
neglects vertical gradients and the simple implementation 
of the Hanle-cffcct restricts the reliable determination 
of the azimuthal angle to regions where the magnetic 
field is strongly inclined to the line of sight. Furthermore 
this method neglects the changes in the polarization 
signal caused by the incomplete Paschen-Back splittin g 



( Socas-Navarro. Trujillo Bueno. fc Landi Degl'Innocentil . 
2004). From a preliminary comparison of the obtained 
Ha slit jaw images and the inferred magnetic loops we 
are confident that the retrieved magnetic field topology 
is close to the real situation. In the following we use 
the term "observed loops" to name the loops inferred 
under these assumptions and to distinguish them from 
loops computed with the help of extrapolations from a 
photosphcric magnctogram. 

3. Computation of 3D magnetic fields from 
photospheric magnetic field measurements. 

A number of authors have modelled the coronal magnetic 
field by extrapolating from photospheric magnetic field 
observations. It is generally assumed that the magnetic 
pressure in the corona is much higher than the plasma 
pressure (small plasma (3) and that therefore the magnetic 
field is nearly force-free. The extrapolation methods based 
on t his assumpti on include potential field extrapolation 
/e.g.lSemell.ll967h . linear force- free field extrapolation (e.g. 



19771 : ScchafeJ, 119781 : IScchafcrl. 11982 



Chiu and Hiltor 

Semell. Il988f) and nonlinear force-free field extrapolation 
e. g.lSakura1 . Jl98ll : Roumeliotid.ll996t lAmari et al.l.ll997 : 
119991 Iwhe. ~ 
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20041) . Force-free magnetic fields have to obey the equa- 



tions 

(V x B) x B = 0, 
V B = 0. 
which are equivalent to 
(V x B) = qB, 
B-Va = 0. 



(1) 

(2) 

(3) 
(4) 



In general a is a function of space. Taking this into ac- 
count corresponds to the non-linear force- free approach. A 
popular simplification is to choose a = constant in the en- 
tire computational domain, the linear force-free approach. 
The choice a = corresponds to current-free potential 
fields. In this paper we compute potential fields, linear 
force-free and non-linear force-free fields and compare the 
result with the magnetic loops reconstructed from the ob- 
servations. 

3.1. Potential and linear force-free fields. 



We use the method of ISeehaferl (|1978l ) for calculating 
the linear force-free field. The method requires a, line-of- 
sight magnetogram and contains a free scalar parameter a, 



where a = corresponds to potential fields. The method 
gives the components of the magnetic field in terms of a 
Fourier series. The observed magnetogram which covers 
a rectangular region extending from to L x in x and 
to L y in y is artificially extended onto a rectangular re- 
gion covering — L x to L x and — L y to L y by taking an 
antisymmetric mirror image of the original magnetogram 
in the extended region, i.e. B z (—x,y) — —B z (x,y) and 
B z (x, — y) = —B z (x, y). The advantage of taking the anti- 
symmetric extension of the original magnctogram is that 
the extended magnetogram is automatically flux balanced. 
We use a Fast Fourier Transformation (FFT) scheme to 
determine the coefficients of the Fourier series, a has the 
dimension of an inverse length. As a characteristic length 
scale we choos e the harmonic mean L of L x and L y . (See 
Seehaferlll978l for details.) 



3.2. Non-linear force-free fields. 

We solve Eqs. (Ill) and © by means of an optimiz ation 



principle (|Wheatland et al.ll200ol IWiegelmannll2004h : 



L= / w(x,y,z) [ J B- 2 |(VxB)xB| 2 + |V-B| 2 ] d 3 x, (5) 
Jv 

where w(x, y, z) is a weighting function. It is obvious that 
(for w > 0) the force-free Eqs. (jT][2j) arc fulfilled when 
L equals zero. We compute the magnetic field in a box 
with nx = 95, ny = 68 and nz = 40 points. The numer- 
ical method works as follows. As an initial configuration 
we compute a pote ntial magnetic fi eld in the whole box 
with the help of the lSeehaferl ()1978f ) method. As the next 
step we use photospheric vector magnetic field data to 
prescribe the bottom boundary (photosphere) of the com- 
putational box. On the lateral and top boundaries the field 
is chosen from the potential field above. We iterate for the 
magnetic field inside the computational box by minimiz- 
ing Eq. ([5]) . The weighting function w equals 1 everywhere 
in the computational box except in a boundary layer of 10 
points towards the lateral and top boundary of the com- 
putational box, where w decreases smoothly to with a 
cosine function. The boundary layer diminishes the influ- 
ence of the lateral and top bou ndary conditions o nto the 
magnetic field in the box. (See IWiegelmannl 2004 for de- 
tails.) 

4. Results 

To compare the reconstructed magnetic field with the ob- 
served magnetic loops wc compute magnetic field lines 
from the reconstructed fields using a fourth order Rungc- 
Kutta field-line tracer. The field-line tracer starts the in- 
tegration at any arbitrary point in space and traces the 
magnetic field in the +B and —B direction until the pho- 
tosphere is reached in both directions. As a measure of 
how well the magnetic field lines and the observed loops 
agree, we compute the spatial distance of the two curves 
in 3D integrated along the whole loop length I from r = 
to r = I. As a result we get a dimensionless number 




Bz /G 
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Fig. 1. Observed and extrapolated fields. The left panels show the field lines in 3D and the right panels projections on 
the photosphere. The photosphcric LOS field is represented by colour-coding. Top panels: Original observed loops. We 
show 14 out of a total of 39 loops studied here. Second row: potential field extrapolation. Third row: linear force-free 
field extrapolation with a ■ L = 3.0. Fourth: non-linear force-free field extrapolation. The fifth row shows one loop 
(Black: original, white: potential field, orange: linear force-free field with a-L = 3.0, yellow: non-linear force-free field.) 
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C = ji f \J (r obs (r) - r oxt rapoi (r)) 2 dr, where r is the geo- 
metrical length measured along the loop and C = if both 
curves coinci de. This comparison method has been previ- 
ously used by Wiegelmann fc Neukirchl ( 2002f ) to compare 
magnetic field lines with stereoscopically observed loops. 

We start the field line integration at 20 points each 
chosen to lie along an observed field line r b s (i~) and com- 
pute 20 corresponding field lines. We compare the shapes 
of these computed field lines (loops) with the observed 
field line and compute the quantitative measure C . The 
lowest value of C corresponds to the optimal computed 
magnetic field line. For potential fields and non-linear 
force-free fields the choice of the starting point is the only 
free parameter and finding the optimal field line is a one- 
dimensional minimization problem. Linear force- free fields 
have the free parameter a and computing the optimal lin- 
ear force-free field line is a two-dimensional minimization 
problem with respect to a and the integration starting 
point. 

Fourteen representative loops are shown in Fig. [1] with 
observed, potential, linear and non-linear force-free loops 
being plotted from top to bottom. This figure shows, that 
force-free extrapolations are superior to potential fields, 
with non-linear force-free being apparently the closest to 
the observations. A more quantitative measure is provided 
by the C values, which are given for all 39 observed loops 
in Fig. [5] and Table [TJ The lower the value of C the better 
the reconstructed loops agree with the observed loops. A 
value of C < 0.35 seems to be acceptable. 

We find that the simplest magnetic field model, po- 
tential fields (second row in Fig. [1] Fig. [3] a), provides 
no agreement with any observed loop (top row in Fig. [T]) 
within the C < 0.35 limit. 

The inclusion of electric currents, in lowest order with 
the linear force- free approximation (Fig.[^]b) provides bet- 
ter results, and for 35% of the observed loops we get sat- 
isfying agreement with the observations. One has to keep 
in mind, however, that a consistent linear force-free recon- 
struction requires a unique value of a in the entire con- 
sidered volume. Most of the loops have an optimum value 
of aL in the range 3 — 4 and the quality criterion C only 
changes slightly within this range, so that a unique value 
of a in this range does not give significantly worse results 
than the optimal value of a. The linear force-free fields in 
Fig. [1] third row have been computed with aL = 3. 

The most involved model used here, the non-linear 
force-free approach (fourth row in Fig. [T] and Fig. [5] c), 
gives even better results than the linear force free ap- 
proach. We get a suitable agreement with the observed 
loops for 64% of the loops within the C < 0.35 limit. 
Let us remark that all observed loops which cannot be 
reconstructed with this model have at least one foot point 
close to the boundary of the available vector magnetic field 
data, (see Fig. [3]). 
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Fig. 2. Correspondence of extrapolated loops and mea- 
sured loops for all 39 observed loops. The panels (a)-(c) 
show distributions over C for potential, linear and non- 
linear force-free fields, respectively. Panel (d) shows the 
value for all loops. The rhombi correspond to a potential 
field reconstruction, the stars to the optimal linear force- 
free reconstruction and the triangles to a non-linear force- 
free reconstruction. (See also Table [1] for the C- values of 
the individual loops.) 



5. Conclusions 

We have compared the observationally inferred structure 
of magnetic loops in the upper chromosphere with mag- 
netic fields extrapolated from photosphcric measurements. 
We find that the simplest model, potential fields, is not 
sufficient to reproduce the observations. The inclusion of 
field line parallel currents, the so called force- free ap- 
proach, gives much better results. Among force-free mod- 
els, a linear one gives less accurate results than the non- 
linear force-free extrapolation. 

We find that the observed and extrapolated loops agree 
quite well for almost 2/3 of the loops, while the remain- 
ing 1/3 might suffer from the limited field of view of the 
available vector magnetogram. 

The investigated active region is quite young. With a 
vertical upflow speed of v = 1.5_t0.5km/s at the loop apex 
and a maximum loop height of lOMm we can estimate the 
time elapsed since the loop tops first emerged as 2hT40min 
assuming a constant rise speed. A horizontal shear flow of 
lkm/s on the photosphere would give a shear of 7.2Mm. 
This value is comparable to the difference of the footpoint 
locations between potential field loops and observed loops. 
It is therefore not clear, whether the electric current has 
been caused by shear flow motion or if the magnetic loops 
already contain the current during their emergence. The 
rise of the loops may also explain some of the discrep- 
ancy between observed and extrapolated loops, since the 
loops may change somewhat during the time needed for 
the instrument to scan the region. 
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Table 1. The first column corresponds to an arbitrary loop 
number of the observed loops, the second column compares the 
measured loops with a potential field extrapolation. Cols. 3 and 
4 compare the observed loops with linear force-free fields and 
aL = 3 and aL = 4, respectively. The fifth column contains the 
optimal linear force-free parameter aL and the sixth column 
the comparison with the observed loops for this optimal value. 
In the last column we compare the observations with non-linear 
force-free reconstructed magnetic loops. 



Nr 


C'pot 


Clin 

aL = 3 


Clin 

aL = 4 


«opt 


Clin 

a op t 


C*nonlin 





1.69 


1.40 


1.06 


4.2 


1.04 


1.60 


1 


1.13 


0.83 


0.86 


2.6 


0.79 


0.86 


2 


0.48 


0.30 


0.28 


3.5 


0.27 


0.15 


3 


0.57 


0.34 


0.31 


3.7 


0.30 


0.14 


4 


0.75 


0.64 


0.65 


3.2 


0.63 


0.65 


5 


0.80 


0.52 


0.45 


3.8 


0.45 


0.19 


6 


0.69 


0.37 


0.31 


4.1 


0.30 


0.17 


7 


0.90 


0.56 


0.43 


4.3 


0.40 


0.12 


8 


0.58 


0.35 


0.30 


3.8 


0.29 


0.14 


9 


1.90 


1.57 


1.39 


4.2 


1.05 


1.38 


10 


0.68 


0.44 


0.51 


3.0 


0.44 


0.59 


11 


0.43 


0.26 


0.29 


2.9 


0.25 


0.25 


12 


0.98 


0.56 


0.73 


3.5 


0.52 


1.01 


13 


0.89 


0.86 


0.87 


3.5 


0.85 


0.95 


14 


1.12 


0.66 


0.53 


3.7 


0.52 


0.26 


15 


1.05 


0.63 


0.49 


4.0 


0.49 


0.20 


16 


0.66 


0.35 


0.29 


4.1 


0.29 


0.15 


17 


0.43 


0.27 


0.28 


3.0 


0.27 


0.21 


18 


0.41 


0.26 


0.26 


3.3 


0.25 


0.16 


19 


0.48 


0.30 


0.30 


3.3 


0.28 


0.16 


20 


0.67 


0.35 


0.31 


4.4 


0.29 


0.19 


21 


0.49 


0.34 


0.30 


3.6 


0.29 


0.20 


22 


0.74 


0.53 


0.42 


4.4 


0.42 


0.18 


23 


0.45 


0.31 


0.28 


3.3 


0.27 


0.34 


24 


0.48 


0.35 


0.31 


3.4 


0.31 


0.29 


25 


0.52 


0.38 


0.36 


4.2 


0.34 


0.20 


26 


0.52 


0.38 


0.34 


4.0 


0.34 


0.26 


27 


0.58 


0.40 


0.48 


2.6 


0.37 


0.24 


28 


0.66 


0.45 


0.57 


2.5 


0.39 


0.23 


29 


0.81 


0.61 


0.84 


2.8 


0.60 


0.82 


30 


1.01 


0.68 


1.25 


3.2 


0.67 


1.40 


31 


1.44 


0.98 


0.80 


3.6 


0.63 


2.73 


32 


1.33 


0.70 


0.79 


3.0 


0.70 


0.72 


33 


1.30 


0.71 


0.71 


3.2 


0.64 


0.53 


34 


1.25 


0.72 


0.60 


3.4 


0.57 


0.34 


35 


1.02 


0.64 


0.49 


4.0 


0.49 


0.20 


36 


1.94 


1.66 


1.32 


4.2 


1.30 


1.57 


37 


0.79 


0.48 


0.37 


4.0 


0.37 


0.28 


38 


1.11 


1.27 


1.30 


-0.2 


1.11 


1.50 
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